%mu=0.8 Upwind
clear
diatah=17*5/4.0;
a=(15-diatah):0.05:(25+diatah);
u=f0(a);
n=(10+2*diatah)/0.05+1;
for i=1:25*17
    n=n-2;
    u_temp=zeros(1,n);
    for j=1:n
       u_temp(j)= 0.2*u(j+1)+0.8*u(j);
    end
    u=u_temp;
end
u_temp=u;
figure(1);
axis([15 25 -0.4 1.2]);
set(gca,'XTick',15:1:25);   
set(gca,'yTick',-0.4:0.2:1.2); 

hold on;
k=15:0.05:(25);
z=f17(k);
fp0=plot(k,z);
fp0.Color = 'b';
grid on;
hold on;
fp1=plot(k,u_temp);
fp1.LineStyle='-';
fp1.Color='r';
fp1.Marker='.';
string='Upwind-mu=0.8.png';
title(string);
saveas(1,string);
clf(1);


